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Abstract 

We discuss the evidence for the presence of QCD saturation effects in the data 
collected in d+Au collisions at RHIC. In particular we focus our analysis on 
forward hadron yields and azimuthal correlations. Approaches alternative to 
the CGC description of these two observables are discussed in parallel. 
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1. Introduction 

The large amount of experimental data collected at RHIC during the last 
decade over a wide kinematical range allows to explore novel QCD effects. 
Indeed RHIC measurements present a number of features which are well de- 
scribed within the Color Glass Condensate effective theory of QCD at high 
energies (see e.g. [H, 0] and references therein). The main physical ingredient 
in the CGC is the inclusion of unitarity effects through the proper consider- 
ation of both non-linear recombination effects in the quantum evolution of 
hadronic wave functions and multiple scatterings at the level of particle pro- 
duction. Such effects are expected to be relevant when partons with a small 
enough energy fraction x are probed in the nuclear (or, in full generality, 
hadronic) wave function. 



Email addresses: javier.lopez-albacete@cea.fr (Javier L. Albacete), 
cyrille .marquetOcern. ch (Cyrille Marquet) 



Preprint submitted to Nuclear Physics A 



September 17, 2010 



The CGC formalism itself is only reliable at small x, in that kinematic 
regime gluon densities are large and gluon self-interactions become highly 
probable, thus taming, or saturating, further growth of the gluon occupation 
numbers. While the need for unitarity effects comprised in the CGC is, at a 
theoretical level, clear, the real challenge from a phenomenological point of 
view is to assess to what extent they are present in available data. Such is 
a difficult task, since different physical mechanisms concur in data, and also 
because the limit of asymptotically high energy in which the CGC formalism 
is developed may not be realized in current experiments. 

Notwithstanding these difficulties, we argue below that RHIC data [3|,|4],|5[ 
offer compelling evidence for the presence of saturation effects. Such claim 
is based on the successful simultaneous description of the suppression of 
particle production j^j and azimuthal correlations [?} at forward rapidities 
in d+Au collisions compared to p+p collisions, using the most up-to-date 
theoretical tools available in the CGC approach. We focus on forward parti- 
cle production for the following reason: such processes are sensitive only to 
high-momentum partons inside one of the colliding hadron, which appears 
dilute, while mainly small-momentum partons inside the other dense hadron 
contribute to the scattering. Since the high— x part of a proton wave func- 
tion is well understood in perturbative QCD, forward particle production in 
high-energy d+Au (or generically p+A) collisions is ideal to investigate the 
small— x part of the nucleus wave function. 

In the case of single-inclusive hadron production pA^-hX, denoting p± 
and y the transverse momentum and rapidity of the final state particles, the 
partons that can contribute to the cross section have a fraction of longitudinal 
momentum bounded from below, by x p (for partons from the proton wave 
function) and xa (for partons from the nucleus wave function) given by 

x p = xp , xa = xp e~ 2y , (1) 

where we introduced the Feynman variable xp = \p±\e y / s/snn with ^/snn 
the collision energy per nucleon. With ^/snn 3> \p± I and forward rapidities 
y > 0, the process features x p < 1 and X4<1, meaning that the scattering in- 
volves on the proton side, quantum fluctuations well understood in QCD, and 
on the nucleus side, quantum fluctuations whose non-linear QCD dynamics 
can be studied. 

In the case of double-inclusive hadron production pA-^-h 1 h 2 X, denoting 
Pi±, P2± and yx, y% the transverse momenta and rapidities of the final-state 
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particles, the Feynman variables are Xi = \pi±\e yi /y/s^N an d x p an d x a read 

x p = X\ + x 2 , x A = X\ e~ 2yi + x 2 e~ 2y2 . (2) 

Therefore, only the forward-forward case is sensitive to values of x as small as 
in the single particle production < 1 and xa «C 1 . The central-forward 

measurement does not probe such kinematics: moving one particle forward 
increases significantly the value of x p compared to the central-central case 
(for which x v = x a = \pi\l \J snn), but decreases xa only marginally For 
this reason we shall refer to these two situations as mid-rapidity ones. The 
approximations made in CGC calculations apply best with forward-rapidity 
observables, at RHIC mid-rapidity ones are contaminated too much by large- 
x physics to be treated at a quantitative level. 

Saturation-based approaches were the only ones to correctly predict the 
suppression of forward particle production jsl, 0] and the azimuthal decorre- 



lation of forward hadron pairs [10J. In the following, we also comment on 
alternative mechanisms that successfully describe mid-rapidity particle pro- 
duction in d+Au collisions. We note that, while some of these approaches are 
also able to describe, a posteriori, the suppression of forward yields, we are 
not aware of any formalism that can also describe the azimuthal decorrelation 
of forward hadron pairs. 



2. Nuclear modifications at mid-rapidity 

In relativistic heavy-ion collisions, nuclear effects on single particle pro- 
duction are typically evaluated in terms of ratios of particle yields called 
nuclear modification factors: 

h dN^ hX /dyd 2 p ± 

pA N coll dNPP^x /dyd 2 p± > I ) 

where N co u is the number of nucleon-nucleon collisions in the p+A colli- 
sion. If high-energy nuclear reactions were a mere incoherent superposition 
of nucleon-nucleon collisions, then the observed R p a would be equal to unity. 
However, RHIC measurements in d+Au collisions (or peripheral Au+Au col- 
lisions) [sl, 0] exhibit two opposite regimes: at mid rapidities the nuclear 
modification factors exhibit an enhancement in particle production at in- 
termediate momenta \p±\ ~ 2 4- 4 GeV. In turn, a suppression at smaller 
momenta is observed. However, at forward rapidities the Cronin enhance- 
ment disappears, turning into an almost homogeneous suppression. 
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Figure 1: RHIC experimental results for RJau at 77 = compared with different calcula- 
tions. From left to right: comparison with EPS09 parametrization 
calculation pjj], and the CGC approach of fl3l ] . 



11], a Glauber-Eikonal 



Single-hadron production data at mid-rapidity have been successfully an- 
alyzed through different formalisms and techniques. Below we sketch an 
incomplete, but representative list of the variety in the theory spectrum: 

• Leading-twist perturbation theory. The assumption is that stan- 
dard collinear factorization holds in nuclear reactions, meaning that 
highly- virtual partons in nuclei behave independently as they do in pro- 
tons. For each parton species % the nuclear parton distribution func- 
tions are taken to be proportional to that of a proton: f^(x,Q 2 ) = 
Rf(x, Q 2 ) fi^ix, Q 2 ). The proportionality factors Rf(x, Q 2 ) are fitted, 
in part, to available d+Au data and in some cases, as in the EPS 
parametrization [ll|, evolved according to DGLAP evolution. The re- 
sulting data description is displayed in Fig. [Ik. 

• Glauber-eikonal multiple scatterings. This approach takes into 
account power corrections to the leading-twist approximation. It re- 
lies on a resummation of incoherent multiple scatterings. Typically, 
this results in a momentum broadening of the scattered parton which 
is responsible for the Cronin enhancement and, due to unitarity con- 
straints, to a depletion of particle production at small transverse mo- 
menta, in agreement with the qualitative features of the data as can 
be seen in Fig. [Tb. Performing the complete resummation including 
energy-momentum conservation is a challenging task. Sometimes, a 
detour of the strict calculation is taken by resorting to unintegrated par- 
ton distributions which include information about the intrinsic trans- 
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verse momentum of the partons ko, and the average transverse momen- 
tum gained during the interaction with the nucleus Ak: f^(x, Q 2 ) — > 
F t A (x,Q 2 ,< kl > + < Ak 2 > (y/s,b,p±). While the intrinsic k is 
adjusted in p+p collisions, the gained transverse momentum is let to 
depend on the collision energy, centrality and p± of the detected hadron. 

• Color Glass Condensate. The CGC approach resums all power cor- 
rections which are dominant in the high-energy/small-x limit. It relies 
on two main assumptions: the scattering process is fully coherent and 
the propagation of the leading parton through the nucleus is eikonal, i.e. 
the momentum transfered during the collision is only transverse. Then, 
non-linearities or saturation effects can be taken into account either 
at the semiclassical level or, more accurately, including the quantum 
corrections embodied in the JIMWLK evolution equation. The work 
presented in Fig. [lb relies on a combination of both, with quantum 
corrections modeled according to analytical estimates [13(. Note that 
in the CGC framework, the saturation scale Q s which characterizes the 
onset of non-linear effects in the nuclear wave function, also determines 
both the intrinsic transverse momentum and the amount of it gained 
during the collision. In the small-a: limit, it is not possible to distin- 
guish saturation from multiple scatterings, such a separation would be 
frame dependent. Both become important when a large gluon density 
is reached, and including one without the other is not consistent. 

Simply by looking at Fig. 1, one concludes that the three different ap- 
proaches above offer a comparably good description of the data, so it is 
difficult to extract any clean conclusion about the physical origin of the 
Cronin enhancement. This is probably due to the kinematic region probed 
in these measurements. For a hadron momentum of \p±\ ~ 2 GeV, one 
is sensitive to (x p ~ 0.01 -j- 0.1. In this region different physical 

mechanisms concur, so neither of the physical assumptions underlying the 
approaches above is comletely fulfilled. Indeed, the coherence length, esti- 
mated to be l c ~ l/(2mjv xa) ~ 1 -S- 10 fm is, on average, smaller than the 
nuclear length, so the fully coherent treatment of the collision implicit in the 
CGC is not completely justified, and large- a; corrections are expected to be 
relevant. Moreover, both the coherent or incoherent treatments of the colli- 
sion in 12|, |l3| need to invoke the presence of an intrinsic scale, presumably 
of non-perturbative origin, of the order of 1 Gev in order to reproduce the 
data, whose dynamical origin is not evident at all. 
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Figure 2: The coincidence probability at mid-rapidity as a function of A(j>. RHIC data 
show that in the central-central case the away-side peak is similar in d+Au and p+p 
collisions. In the central-forward case, the Glauber-eikonal (left plot from [l3|) and CGC 
(right plot from [lU) calculations predict that the away-side peak is suppressed in d+Au 
compared to p+p collisions, in agreement with later data. 



More insight is gained with di-hadron correlations. In particular, let 
us focus on the A0 dependence of the double-inclusive hadron spectrum, 
where A<f> is the difference between the azimuthal angles of the measured 
particles hi and hi- Nuclear effects on di-hadron correlations are typically 
evaluated in terms of the coincidence probability to, given a trigger particle 
in a certain momentum range, produce an associated particle in another 
momentum range. In a p+p or p+A collision, the coincidence probability is 
given by OP(A0) = N pair (A(f))/N trig with 

f dN pA ^ hlh2X f dN pA ~* hX 

N *°^ = J j„jp, • N <- = J • (4) 

yi,\Pi±\ v, px 

First measurements were performed at RHIC at mid-rapidity by the PHENIX 
and STAR collaborations 0, 16]. In the central-central case, the coincidence 



probability features a near-side peak around A0 = 0, when both measured 
particles belong to the same mini-jet, and an away-side peak around A<p = 
7r, corresponding to hadrons produced back-to-back. In the central-forward 
case, there is naturally no near-side peak. 
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Associated: 3.i<T\< 3.9, pT - 0.45-1. 59 GeV/c 
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Figure 3: Central- forward preliminary IdAu data as a function of centrality 

;i7]. 

In central 

collisions, the integral of the coincidence probabilty is about half that in p+p collisions, 
reflecting the depletion of the away-side peak. Collinear factorization cannot reproduce 
this behavior, while both Glauber-eikonal and CGC calculations predicted it. 



Either in p+p or d+Au collisions, the sizeable width of the away-side peak 
cannot be described within the leading-twist collinear factorization frame- 
work. This indicates that, while it may not be obvious in single particle 
production, power corrections are important when \p±\ ~ 2 GeV. At such 
low transverse momenta, collinear factorization does not provide a global 
picture of particle production at RHIC, even at mid-rapidity. On the con- 



trary, both the Glauber multiple scattering [14| and CGC [15| approaches 
can qualitatively describe the data, including the depletion of the away-side 
peak in d+Au collisions when going from central-central to central-forward 
production. This is illustrated in Fig. [2j Such a depletion does not occur in 
p+p collisions, it is due to nuclear-enhanced power corrections, and therefore 
the p+A to p+p ratio of the integrated coincidence probabilites 

/rf0 C7y(A0) 

is below unity. In Fig. [3j recent PHENIX data on I d A u are displayed as a 
function of centrality. At the moment, since xa is not that small, it is not 
clear whether the mechanism for this suppression is due to saturation effects 
rather than incoherent multiple scatterings. 
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3. Moving forward 



As outlined in the introduction, data collected in the deuteron fragmen- 
tation region offer a cleaner opportunity to explore saturation effects. Let 
us first explain our implementation of the CGC framework. The CGC is 
equipped with a set of non- linear renormalization group equations to describe 
the evolution of hadronic wave functions towards small-x. In the large- N c 



limit, they reduce to the BK equation [18l. Il9f. The recent determination of 
running coupling corrections to the original leading-log equations 20|, |2l|] has 
proven an essential step in promoting the BK equation to a phenomenological 
tool. Indeed, the running coupling BK equation (rcBK) has been employed 



to successfully describe inclusive structure functions in e+p scattering [22 



and also the energy and multiplicity dependence of total hadron multiplicities 
in Au+Au collisions at RHIC 
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The rcBK equation reads 



Af(r, Y) 
(91n(l/x) 



J d 2 ri # run (r, ri) r 2 ) [M(rt,Y)+M{r 2 ,Y) 
-J\f(r, Y) - J\f(n, Y) J\f(r 2 , Y)} , 



(6) 



where r 2 = r — ri (we use the notation v = |v| for two-dimensional vectors 
in Eq. ([5]) and Eq. ([7])). Af(r,Y) is the dipole scattering amplitude in the 
fundamental representation, with Y = \n(xo/x) the rapidity and r the dipole 
transverse size. It turns out that the evolution kernel 



K T 



vr,ri,r 2 J 



N c a s {r 2 ) 
2tt 2 



L r l 



a Art 



-1 + 



1 



(7) 

proposed in |20| minimizes the role of higher order corrections, making it 
better suited for phenomenological applications. Detailed discussions about 
other prescriptions proposed to define the running coupling kernel, and about 
the numerical method to solve the rcBK equation can be found in 24 . 

Eq. <Q needs to be suplemented with initial conditions, which we choose 
to be of the McLerran-Venugopalan type: 



Af(r,Y=0) 



exp 



Ar 



+ e 



(8) 



where A = 0.241 GeV. This introduces two free parameters: the value Xq 
where the evolution starts and the initial saturation scale felt by quarks Q s q. 
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3.1. Nuclear modification factors 



According to Ref. [25|, the differential cross section for forward hadron 
production in p+A collisions is given by 



dN pA^hX 

dy d 2 p± 



K Y1 [ ~T x ±fq/p(xpiP±) F ( x Ai— ) D h/q (z,p 2 ± ) 

q Jx F z L \ Z / 

+ x 1 f g/p (x p ,p 2 L ) F (x A , —J D h/g (z,p 2 ± ) 



(9) 



where the unintegrated gluon distributions F and F are related to the dipole 
scattering amplitude through Fourier transformations: 



F(x, k) 
F(x,k) -- 



d 2 r 
(2^)2 
d 2 r 

WT 2 ' 



-ik-r 



-ikr 



[1 -M(r,Y = \n(x /x))] 
[1 -N(r,Y = \n(x /x))} 2 



(10) 
(11) 



In principle F is the Fourier transform of the dipole scattering amplitude 
in the adjoint representation 26|, |27J, l28j . we have used to large- N c limit in 
Eq. ( TTTj) . In Eq. ([9l), fi/ p and Dh/i refer to the parton distribution function 
(pdf) of the incoming proton and to the final-state hadron fragmentation 
function respectively. Here we will use the CTEQ6 NLO pdf 's 29] and the 
DSS NLO fragmentation functions 3(], 31 ]. We have assumed that the factor- 
ization and fragmentation scales are both equal to the transverse momentum 
of the produced hadron. Note that the projectile in our calculations is actu- 
ally a deuteron, and we obtain the neutron pdf from the proton one assuming 
isospin symmetry. 

For light hadron production discussed here, the difference between the 
rapidity and pseudo-rapidity r\ of the produced hadron can be neglected, 
yielding the following kinematics: x p = Xf/z and xa = {xp/z) exp (— 2y) 
with x F = \Jm 2 h + p\j 'y/s NN exp (77) fa \p±\/sJs NN exp (y) introduced be- 
fore. Due to parton fragmentation, the values of x's actually probed are 
generically higher than x p and xa defined in Eq. ([T]). 

With this set up we reach a very good description [6| of the forward 
negatively charged hadron and neutral pion yields measured by the BRAHMS 
[3| and STAR collaborations respectively, in p+p and minimum bias d+Au 
collisions. The parameters of the fit are Xq = 0.025 -j- 0.005 (0.0075) and 
Q 2 = 0.5 -r 0.4 (0.2) GeV 2 for the initial nucleus (proton) wave function. 
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Figure 4: Negatively charged hadrons and neutral pions at forward rapidities at RHIC in 
p+p (left) and minimum bias d+Au (right) collisions, compared to our calculation Q. 

Moreover, no .fT-factors are needed except for the most forward pions rj > 3, 
where we find that K = 0.3 (0.4) is needed to describe the data. 

Our results are shown in Fig. HI By simply taking the ratios of the corre- 
sponding spectra, we get a very good description of the nuclear modification 
factors at forward rapidities, and this is shown in Fig. [5h. It should be noted 
that we use the same normalization as experimentalists do in their analysis 
of minimum bias d+Au collisions: N co u = 7.2. Physically, the observed sup- 
pression is due to the relative enhancement of non-linear terms in the small- x 
evolution of the nuclear wave function with respect to that of a proton. 

After the data confirmed the CGC expectations, it has been argued that 
the observed suppression of particle production at forward rapidities is not 
an effect associated to the small values of xa probed in the nuclear wave 
function, but rather to energy-momentum conservation corrections relevant 



for xf — > 1 |32|, |33j. Such corrections are not present in the CGC, built upon 
the eikonal approximation (this may explain why a X-factor is needed to 
describe the suppression of very forward pions). The energy degradation of 
the projectile parton, neglected in the CGC, through either elastic scattering 
or induced gluon brehmstralung would be larger in a nucleus than in proton 
on account of the stronger color fields of the former, resulting in the relative 
suppression observed in data. A successful description of forward ratios based 



on the energy loss calculation in [32j is shown in Fig. [5b. Calculating forward 
di-hadron correlations in both frameworks could help pin down which is the 
correct picture. In the following, we show that the CGC calculations predicts 
correctly the azimuthal distribution. 
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Figure 5: Nuclear modification factors at forward rapidities in minimum bias d+Au colli- 
sions in CGC [f| (left) and energy-loss 32 1 (right) calculations. 



3.2. Di-hadron azimuthal correlations 

The kinematic range for forward particle detection at RHIC is such that 
x p ~ 0.4 and xa ~ 10~ 3 . Therefore the dominant partonic subprocess is 
initiated by valence quarks in the proton and, at lowest order in a s , the 
pA — > hxhiX double- inclusive cross-section is obtained from the qA — > qgX 
cross-section, the valence quark density in the proton f q / p , and the appropri- 
ate hadron fragmentation functions Dh/ g and Dh/ g : 



dN pA ^ 



hih 2 X 



1 pi /-.I 

dz\ \ dz2 

XI J X2 



*1 I x 2 



X 



dx f q / P (x,fi 2 ) 

"■2 

dNqA ^ q9 x Lp Pl P^ D hl/g ( Zl ,^)D h2/g (z 2 ,^ 

V z l z 2, 



dNQ A^ qg X I ^ P2 f Pi j Dhi/g ( Zl ^*) Dh2/q ( Z2j ^ 
z 2 z l 



.(12) 



We shall use CTEQ6 NLO quark distributions |29| and KKP NLO fragmen- 



tation functions [34J. The factorization and fragmentation scales are both 
chosen equal to the transverse momentum of the leading hadron, which we 
choose to denote hadron 1, p, = \pi±\. We have assumed independent frag- 
mentation of the two final-state hadrons, therefore Eq. (fl~2l cannot be used to 
compute the near-side peak around A$ = 0. Doing so would require the use 
of poorly-known di-hadron fragmentation functions, rather we will focus on 
the away-side peak around A$ = n where saturation effects are important. 
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For the proton, one has x p < x < 1, and if x p would be smaller (this 
will be the case at the LHC), the gluon initiated processes g A — > qqX and 
gA — > ggX should also be included in Eq. f JT2|) . For the nucleus, we shall 
see that the parton momentum fraction varies between xa and e~ 2yi + e~ 2y2 . 
Therefore with large enough rapidities, only the small— x part of the nuclear 
wave function is relevant when calculating the qA — > qgX cross section, and 
that cross section cannot be factorized further: dN qA ^ qgX ^ f g /A®dN q9 ^ q9X . 
Indeed, when probing the saturation regime, dN qA ^ q9X is expected to be a 
non-linear function of the nuclear gluon distribution, which is itself, through 
evolution, a non-linear function of the gluon distribution at higher x. 

Using the CGC approach to describe the small— x part of the nucleus 



wave function, the qA — > qgX cross section was calculated in [35|, |36|, [37, 



10j. It was found that the nucleus cannot be described by only the single- 
gluon distribution, a direct consequence of the fact that small-x gluons in the 
nuclear wave function behave coherently, and not individually. The qA — > 
qgX cross section is instead expressed in terms of correlators of Wilson lines 
(which account for multiple scatterings), with up to a six-point correlator 
averaged over the CGC wave function. At the moment, it is not known 
how to practically evaluate the six-point function. In the large- iV c limit, it 
factorizes into a dipole scattering amplitude times a trace of four Wilson 
lines, and the latter can in principle be obtained by solving an evolution 



equation written down in [35] • However, this implies a significant amount of 
numerical work and also requires to introduce an unknown initial condition. 

Rather we shall follow the approach of [lOj, and use an approximation that 
allows to express the six-point function in terms of the two-point function 
Eq. fflOl) . The resulting cross section for the inclusive production of the 
quark-gluon system in the scattering of a quark with momentum xP + off the 
nucleus A reads llOfl : 



dN qA^qgX 



F 



5(xP + -k + -q + ) F(x A ,A) 



d 3 kd 3 q An 2 

X^2\l^(z,k ± -A;x A )-i/;^(z,ks.-zA)\ 2 , (13) 

Xa/3 

where q and k are the momenta the quark and gluon respectively, and with 
A = k± + q± and z = k + /xP + . In this formula, xa denotes the longitudinal 
momentum fraction of the gluon in the nucleus, and xa = X\ e~ 2yi / ' Z\ + 
xi e~ 2y2 /z2 > xa when the cross section Eq. ( TT31) is plugged into Eq. ( 1121) . 
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Figure 6: The A<j> distribution (1 / cr)da/dA(fi computed in at fixed y\ — 3.5, yi = 2 
and \pi±\ = 5 GeV. The azimuthal correlation is suppressed as \p2±\ decreases, this is due 
to the onset of non-linearities in the nuclear wave function. 

The second line of Eq. ( fl3l) features a fc^-factorization breaking term: 

I^(z,k ± ;x) = J d 2 q ± ^(z,q ± )F(x,k ± -q ± ) , (14) 

and where i/j^b * s the well-known amplitude for q —¥ qg splitting (A, a and 
(3 are polarization and helicity indices). While no additionnal information 
than the two-point function is needed to compute Eq. f ll3jl . since higher-point 
correlators needed in principle have been expressed in terms of F(x, k±), the 
cross section is still a non-linear function of that gluon distribution, inval- 
idating k?— factorization. The rather simple form of the k?— factorization 
breaking term is due to the use of a Gaussian CGC color source distribution, 
and to the large— N c limit. Finally, the factor 5(xP + — k + — q + ) in Eq. ( TT3l) 
is due to the fact that the eikonal approximation, valid at high energies, is 
used to compute the qA — > qgX cross section. The delta function is a mani- 
festation of the fact that in a high-energy hadronic collision, the momentum 
transfer is mainly transverse. 

Before comparing our predictions to the recent RHIC data, we first display 
in Fig. |6]the predictions for the A0 distributions made in [lOj when only the 
leading-order BK evolution was known. The different curves are obtained 
for fixed y 1 , y 2 and \pi±\, while \p2±\ is varied. Although these results are 
only qualitative since for instance parton fragmentation was not included, 
the suppression of the A0 distribution as \p%±\ gets closer to the saturation 
scale was predicted. 
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Figure 7: The coincidence probability as a function of for p+p (left) and central d+Au 
(right) collisions. These preliminary data from [j| show a striking nuclear modification of 
di-hadron azimuthal correlations. The away-side peak, corresponding to hadrons emitted 
back-to-back, is prominent in p+p collisions but is absent in the central d+Au case. Such 
production of monojets was anticipated in the CGC as a signal of parton saturation. 



Then, the derivation of the rcBK equation made robust quantitative cal- 
culations possible. Indeed, after the parameters of the theory (x and Q s0 
for both the proton and the gold nucleus) were contrained by single-inclusive 
forward spectra, parameter-free predictions for the coincidence probability 
CP(A(f)) could be made Q. In Eq. (TJJ, when computing N pair from Eq. ( fl"2l) 
and N trig from Eq. (jSJ) for d+Au and p+p collisions, the integration bounds 
for the rapidities are 2.4 < y < 4, in order to compare with RHIC data 
(^J. In addition, such a restriction does insure that only small-momentum 
partons are relevant in the nucleus wave function, as is assumed in our calcu- 
lation. For the trigger (leading) particle \pi± \ > 2 GeV and for the associated 
(sub- leading) hadron 1 GeV < \p2±\ < \Pi±\- 

The recent data obtained by the STAR collaboration for the coincidence 
probability obtained with two neutral pions are displayed in Fig. [7] for both 
p+p and central d+Au collisions. The nuclear modification of the di-pion 
azimuthal correlation is quite impressive, the prominent away-side peak in 
p+p collisions is absent in central d+Au collisions, in agreement with the 
behavior predicted in Il0|. In Fig. [SK, these data are compared with the 
rcBK calculations of [7|. As mentioned before, the complete near-side peak 
is not computed, as Eq. (fT2l does not apply around A<fi = 0. 
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Figure 8: The coincidence probability as a function of A<f). Left: for p+p and central 
d+Au collisions, preliminary data are compared with CGC predictions; the away-side 
peak in p+p is qualitatively described by the CGC calculation while its disapearence in 
central d+Au is quantitatively consistent with the prediction. Right: CGC predictions 
for different centralities of the d+Au collision; the near-side peak is independent of the 
centrality, while the away-side peak reappears as collisions are more and more peripheral. 

We see that the disapearence of the away-side peak in central d+Au colli- 
sions, compared to p+p collisions, is quantitatively consistent with the CGC 
calculations. The latter are only robust for the d+Au case, but the extrap- 
olation to the p+p case is displayed in order to show that it is qualitatively 
consistent with the presence of the away-side peak in p+p, and also with the 
fact that the near-side peak is identical in the two cases, and is not sensitive 
to saturation physics. Note that since uncorrelated background has not been 
extracted from the data, the overall normalization of the data points has 
been adjusted by subtracting a constant shift, as indicated on the figure. 

To deal with the centrality dependence, we have identified the centrality 
averaged initial saturation scale C) 2 , extracted from minimum-bias single- 
inclusive hadron production data, with the value of Q 2 s0 at b = 5.47 fm, and 
used the Woods-Saxon distribution TaQ>) to calculate the saturation scale at 
other centralities: 

Q 2 M = ^j? , Q 2 = 0.4GeV 2 . (15) 

The value used in Fig. in the central d+Au case is C} 2 (0) ~ 0.6 GeV 2 
at xo = 0.02. The corresponding saturation scale felt by gluons is about 1.2 
GeV 2 and of course it gets bigger with decreasing x. 
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In Fig. [Sb, we show the centrality dependence of the coincidence proba- 
bility. Although it is difficult to trust our formalism all the way to peripheral 
collisions, we predict that the near-side peak does not change with centrality, 
and that the away-side peak reappears for less central collisions. This is con- 
sistent with the fact that peripheral d+Au collisions are p+p collisions. The 
fact that the away-side peak disappears from peripheral to central collisions 
shows that indeed this effect is correlated with the nuclear density. 

Moreover dihadron correlations at mid-rapidity, which are sensitive to 
larger values of xa, feature an away-side peak whatever the centrality. The 
fact that for central collisions the away-side peak disappears from central 
to forward rapidities also indicates that the effect is correlated with the nu- 
clear gluon density. In a similar way, we predict that for higher transverse 
momenta, the away-side peak will reappear. 

We are not aware of any descriptions of this phenomena that does not 
invoke saturation effects. We note that apart from our CGC calculation, a 
successful description based on the KLN saturation model was also recently 



proposed [38]. There, although different assumptions are used, the existence 
of the saturation scale is the crucial ingredient to successfully reproduce the 
data. While more differential measurements of the coincidence probability, 
as a function of transverse momentum or rapidity, will provide further quan- 
titative tests of our CGC predictions, the piece of evidence we have discussed 
in this work strongly indicates that we have observed manifestations of the 
saturation regime of QCD at RHIC. 



4. Predictions for the LHC 

The huge leap forward in collision energy reached at the LHC allows for 
an exploration of small- x effects already at mid-rapidity. There, both the 
target and projectile are probed at small values of x, and energy loss effects 
associated to large-x^ effects are expected to be small. In Fig. Owe present 
our CGC predictions for the nuclear modification factor for negative charged 
hadrons in p+Pb collisions at two LHC energies. Our curves correspond 
to rapidities y = 2 and larger. Technical difficulties related to the intrinsic 
asymmetry of the hybrid formalism used for particle production prevent us 
from calculating the ratios at mid-rapidity. However, the smooth rapidity 
dependence suggests that a large suppression ~ 0.6 is also expected at mid- 
rapidity in the LHC. It should be taken into account that the normalization 
taken to produce the curves in Fig. [9] was N co u = 3.6. 



16 



1 

0.9 
+ h Q.0.8 



DC 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1, 



Vs NN =8.8TeV y=2,4,6 



N coll =3.6 




hybrid formalism 
kt-factorization 



10 15 20 

P T (GeV) 



1 

0.9 

tb Q.0.8 



DC 



25 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1, 



TeV y=2,4,6 



N coll =3.6 




hybrid formalism 
kt-factorization 



10 15 20 

P T (GeV) 



25 



Figure 9: CGC predictions for the nuclear modification factor in p+Pb collisions at two 
LHC energies and rapidities y = 2, 4 and 6. 



We also compare the y = 2 and 4 curves with predictions obtained with 
the /c^-factorization formalism, in order to check the validity of that ap- 
proach, and especially to test up to what value of y it can be used. The 
&T-factorization formula (see 0]) is valid when the dominant contributions 
to the cross section come from small values of x, for both the projectile 
(x p <C 1) and the target (xa <C 1)- For instance, it only includes gluonic 
degrees of freedom. This approach is clearly insufficient at very forward 
rapidities or large p±, where valence quarks of the projectile are important 
(x p — > 1). However, as can be seen in Fig. [HI both formalisms give comparable 
results, as the lines from &/r-factorization overlap with the uncertainty bands 
spanned by the results from the hybrid formalism. This seems to identify a 
kinematical window where both approximations are valid. 
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